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ABSTRACT 

We report results from a long-term monitoring data set obtained with the Rossi X- 
ray Timing Explorer (RXTE) from the bright globular cluster binary source 4U 1820-30. 
During this monitoring, the source intensity varied from a flux of 2.9 x 10~ 9 ergs cm~ 2 s~ x 
to 8.5 x 10~ 9 ergs cmT 2 s~ 1 in the 2-10 keV band, which correspond to, respectively, a 
count rate of 1,600 cps and 4,500 cps in the nominal 2-60 keV RXTE-PCA band. Its 
energy spectral state was observed to encompass both the banana branches and the 
island state. Kilohertz QPOs were observed to exist in the upper-banana and lower- 
banana branches, as well as in the island state. In particular, we report that the kHz 
QPOs were present only when the soft-color (defined as the flux ratio of the 3.9-5.0 keV 
band to the 1.7-3.9 keV band) was below 0.78 without regard to the overall count rate 
and without regard to whether the source is in any one of the three branches or states. 
Their centroid frequency was well correlated with the overall count rate below 2,500 cps. 
When the count rate was above 2,500 cps, the QPO frequency was consistent with a 
constant independent of the overall count rate. We take this as further evidence that 
the inner edge of the accretion had reached the marginally stable orbit, and therefore 
the constant 1060±20 Hz is the marginally stable orbit frequency. 

Subject headings: accretion, accretion disks — stars: neutron — starsdndividual (4U 1820- 
3 ) — X-rays : st ar s 
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1. Introduction 



most important features of our results. 



Low mass X-ray binaries are enigmatic objects. 
They are the brightest X-ray sources in the sky. There 
is little doubt that they harbor fast-spinning neutron 
stars, yet all efforts to date of searching for coherent 
pulsations in their persi stent flux have resulte d in neg- 
ative results (see, e.g. Vaughan et al. 1994 ). Those 
searches, however, have uncovered quasi-periodic flux 
oscillations which are belie ved to be indirectly related 
to the neutron star spin (Alpar & Shaham, van dcr 



Klis et al. 1985 ). Most rece ntly, since their discovery 
with RXTE i n 4U 1728-34 QStrohmayer et al. 19961 ) 
and Sco X-l (van der Klis et al. 1996), the highest 
frequencies observed from any cosmic sources, the so- 
called kilohertz quasi-periodic oscillations, have been 
observed from nearly every LMXB. For a comprehen- 
sive review of the literature on this subject, see van 
der Klis (1997). 

The bright globular cluster binary 4U 1820-30 is 
one of the most studied low mass X-ray binaries be- 
cause of its membership in a globular cluster, its 
Type-I X-ray bursts ( Grindlay et al. 1976] ), its vari- 
ability on all time scales, its short orbital period of 
only 685 s, and most recently, the kilohertz quasi- 
periodic oscillations in the persistent flux ( [Smale ei 
al. 1997). Since it is in the globular cluster NGC 
6624, its distance of 6.4 ± 0.6 kpc is probably more 
accurately known to us than any other LMXB. Its 
short orbital period — the shortest known of any cos- 
mic binary — strongly suggests that its mass-donating 
star is most l ikely a helium dwarf w ith 0.06 to 0.08 
solar masses ( Rappaport et al. 1987 ). 

The characteristics of the kilohertz QPOs d iscov- 
ered in the persistent flux ( Smale et al. 1997 ) were 
rather similar to those of other bursters. Their cen- 
troid frequency was quite tightly correlated with the 
2-60 keV count rate. Their fractional rms amplitudes 
were of the order of 3%. The frequency difference of 
the twin QPO peaks was consistent with a constant 
275 ± 8 Hz. Their fractional rms amplitudes increase 
dramatically with energy; the 8-20 keV amplitudes 
are 1.8-2.9 times greater than those measured in the 
1-8 keV band. 

In this paper we report on data sets that sample 
the large count rate and energy spectral span of this 
source over a calendar period of 6 months starting in 
February and ending in September of 1997. In the 
discussion section, we place our results in the overall 
context of kHz QPO phcnomcnologies and discuss the 



2. Observations 

The RXTE satellite was launched on 30 December 

1995. It carries three instruments: an array of 5 xenon 
gas proportional counters with a nominal energy band 
of 2-60 keV, a high energy X-ray timing experiment 
(HEXTE) with a bandwidth of 15-250 keV, and an 
all sky monitor. In this paper we use only the data 
from the proportional counter array (PCA) because 
of its large area and high sensitivity in detecting kHz 
QPOs. 

The data we analyzed for this paper came from two 
sets of observations. The first set with an exposure 
time of 78 ks occurred between 26 and 31 October 

1996. The second set with a total exposure of 153 
ks accumulated between 9 February and 10 Septem- 
ber 1997. The results from the first data set were 
reported previously in the kHz QPO discovery paper 
(Smale et al. 1997). But we include those measure- 
ments here for completeness. The second set of ob- 
servations were designed to sample the source over its 
super-orbital intensity variation period of about 176 
days ( Priedhorsky fc Terrell 1984 ). The combined 
data set of 231 ks was broken up into 90 segments 
by the RXTE planning process, the south Atlantic 
anomaly, and Earth occult of the source, with the 
shortest contiguous exposure of 128 s and the longest 
3,840 s. The typical segment lasted 3,000 s. The 
data were collected in an event mode with a 128/ii 
(/it = 2~ 20 /is) time resolution and a 64 channel en- 
ergy resolution, in addition to the standard one and 
two data. One Type-I X-ray burst was observed. A 
timing analysis of the burst data revealed no oscilla- 
tions similar to those observed in bursts from other 
sources, such as 4U 1728-34 ( Strohmayer et al. 1996| ). 

A comparison with the long term light curve ob- 
tained with the RXTE- ASM indicates that the second 
data set indeed sampled a 176 day cycle of 4U 1820- 
30 quite adequately. It covered both the low and high 
states of the source, whereas the first data set covered 
the low state. The count rate ranges from a low of 
1,500 cps to 4,500 cps in the nominal 2-60 keV of the 
RXTE-PCA. 

3. Data Analysis and Results 

We Fourier analyzed the data of every 128 s of 
observation with a time resolution of 256/ii. Then we 
obtained the average power spectrum for each of the 
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90 segments. In order to display and visually inspect 
each of the 90 power spectra, we rebinned each one of 
them by a factor of 1024 in frequency so as to have a 
bin size of 8 Hz. A fit to one or two Gaussian profiles 
were performed on those spectra which had visually 
identifiable peaks. 

The energy spectral colors were calculated using 
the standard two data. For each of the 90 segments, 
we obtained its energy spectrum from the standard 
two data set with charged particle and diffuse cos- 
mic X-ray background subtracted using the standard 
RXTE-PCA background estimator. 

Figure |l| shows the X-ray color-color diagram. 
Each of the 90 segments is represented by one data 
point. The filled circles are those segments where kHz 
QPO peaks are visually identified and numerically fit- 
ted to Gaussian profiles. The open circles are those 
segments from which no QPO peaks can be identified 
either visually or numerically. It is clear that with 
only 3 exceptions, the open circles occupy the space 
to the right of soft color 0.78 and the filled circles to 
the left. 

In order to ascertain the nature of the energy spec- 
tral states of those 90 segments of data, we have exam- 
ined their power spectra in the frequency range of 0.02 
to 100 Hz in comparison with the classification scheme 
as defined by Hasinger and van der Klis (1989). We 
conclude that the points in Figure |l| encompass both 
the lower banana and the upper banana branches, as 
well as the island state. In particular we note that 
there is not any systematic correspondence between 
these branches/states and the positions in Figure |l| 
because our data cover a large calendar time span 
and it has been known that the position of the source 
in the color-color diagram displays secular shifts. In 
other words, one can not read off the branch or state 
of the source from Figure |l|. As examples, Figures ^ 
show the power spectra for the three points marked 
in Figure |l|, which correspond to, respectively, the 
upper banana branch, the lower banana branch, and 
the island state. On the one hand, we note that kHz 
QPOs are observed for each of the three segments 
whose low frequency power spectra are shown in Fig- 
ures |^. On the other hand, we also note that there 
are other segments with very similar low frequency 
power spectra as shown in the three panels of Fig- 
ures H which do not show any kHz QPOs. Therefore 
we conclude that the existence of kHz QPOs is not 
directly related to whether the source is in the up- 
per banana, lower banana, or island state. Rather 



most likely, their existence is related to the energy 
spectral state as directly defined by the hardness ra- 
tios as shown in Figure [l]. To demonstrate that the 
separation between the open and filled circles in Fig- 
ure [j] is rather sharp, we show in Figure I the FFT 
power spectra that correspond to the two points as 
indicated by the two arrows in Figure Panel (a) 
is from a 10 September 1997 observation with an ex- 
posure time of 1,843 s and a count rate of 3,513 cps 
and (b) from a 9 February 1997 observation with an 
exposure time of 2,632 s and a count rate of 2,979 cps. 
The measured fractional rms amplitudes from the two 
peaks in panel (b) are 3.2%± 0.7% and 2.9%±0.7%, 
respectively. Given the reasonable assumption that 
any QPO in panel (a) should have the same charac- 
teristics as those in panel (b), we set a 90% confidence 
level upper limit of 1.3% which is typical for all other 
segments of data from which no kHz QPOs have been 
positive detected. 

Figures ^ and || show the QPO characteristics 
as functions of the overall source count rate in the 
RXTE-PCA detectors. The centroid frequency, as 
shown in Figures || is well correlated with the count 
rate below 2,500 cps. Above 2,500 cps, The QPO 
frequencies are consistent with constant 1060±20 Hz 
independent of count rate. 

In the entire count rate span of 1,600 to 3,200 cps, 
the lower QPO and the higher QPO centroid frequen- 
cies change in unison and maintain a difference con- 
sistent with a constant as shown in Figure |a. The 
fractional rms amplitude is also correlated with the 
count rate over the entire count rate span as shown 
in Figure ||b. We note here that there is ambiguity 
in determining whether a QPO peak corresponds to 
the lower or the higher QPO of the pair when it is 
the only peak observed. In Figures ^ and 0, we fol- 
lowed the most natural identification in that the open 
circles and the filled circles form more or less smooth 
lines when connected. 

When the count rate increases from 1,600 cps to 
3,200 cps and above, the source state transitions from 
the island state to the banana branches. At the 
counte rate 2,500 cps where the slope of frequency 
vs. count rate changes, we do not observe any discon- 
tinuous change in the 0.2-100 Hz parts of the power 
spectra. The points on the plateau in Figure || are on 
the lower or upper banana branch. 
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4. Discussion 

We have observed kHz QPOs from 4U 1820-30 on 
both the banana branches and the island state. This 
is the first source from which kHz QP Os have been 
ob served in a ll branches or states ( see Wijnands et 
al. If 997b and Wijnands et al. 1997a for discussions). 
The existence of kHz QPOs in the flux is strongly 
correlated with the source position on the color-color 
diagram. We have examined the correlation between 
the source states as defined by Hasinger and van der 
Klis (1989) and the existence of kHz QPOs and found 
that there does not seem to be a one to one correspon- 
dence between the two. For example, we found kHz 
QPOs in all three branches/state in some segments 
of the data and found no kHz QPOs in any of the 
three branches/state in other segments of the data. 
Therefore we think that over long terms color-color 
daigram alone may be a better predictor or indicator 
for the existence of kHz QPOs. 

Figure ^ demonstrates the well-defined correlation 
between the QPO centroid frequency and the overall 
count rate. It is similar to nearly all other sources for 
which this correlation is known. Although the data 
analyzed here covers a large span of calendar time, 
there were no large shifts in this correlation, in con- 



trast to the c ases of 4U 0614 +091 ( |Ford et al 1997a|) 
4U 1608 -52 ( |Yu et al. 1997| ), and Aql X-l QZhang et 
al. 199£ ) . A very important feature of Figure [|, com- 
pared to the frequency vs. count rate correlations 
of other sources, is that the QPO frequency is inde- 



pe ndent of the count rate in the range of 2 , 600 to 



3,200 cps. This is a feature that has been proposed 



(Miller, Lamb, & Psaltis 1998) to strongly support 


the hypothesis (Kaaret, Ford, & Chen 1997, 


Zhang, 


Stohmayer, & Swank 1997 


) that the highest observed 



QPO frequencies are the marginally stable orbit fre- 
quencies near the neutron stars in the low mass X- 
ray binaries. Figure || strongly suggests that indeed 
we have observed the highest QPO frequency from 
4U 1820-30. If this frequency is the marginally sta- 
ble orbit frequency, then neutron star in 4U 1820-30 
would have a mass in the vicinity of 2.2Af Q . Given 
the importance of this conclusion, more corroborating 
evidence from other sources is needed. 

Figure [|b tests the constancy of the frequency dif- 
ference over a much larger span of count rate than 
previously possible flSmale et al. 1997 ) . With the sta- 
tistical accuracy we can measure, the difference is a 
constant 275±8 Hz. This is similar to the cases of 



4U 1728-34 and 4U 0614+091, where the differences 
are consistent with being constants 363 and 327 Hz 
( gtrohmayer et al. 1996] , Ford et al. 1997h| ), respec- 
tively; in contrast to the cases of Sco X-l and 4U 
1608-52 whose frequency differences have been ob- 



served to change by over 100 Hz (van der Klis et al 
1997b|, [Mendez et al. 1997|). 



The fractional rms amplitude as a function count 
rate, Figure |^c, indicates that, on average, when the 
count rate decreases by a factor of 2 from 3,200 cps 
to 1,600 cps, the rms amplitude increases by a fac- 
tor of 2 from about 3% to about 6%. This indicates 
that the flux pulsed at the QPO frequency has a con- 
stant amplitude. It is only the un-pulsed flux that has 
changed by a factor of 2. This point can be poten- 
tially constraining as it requires that the overall X-ray 
luminosity of the source change by as much as a fac- 
tor of 2, but the flux pulsed at the kHz frequencies 
remain more or less the same. 

As the source count rate decreses, the upper fre- 
quency QPO peak becomes broader and its amplitude 
larger. The lower frequency QPO would be sand- 
wiched in between the low frequency noise which is 
below 100 Hz and the upper QPO if it were below 
the upper frequency by 275 Hz. We have no positive 
detection of it and set a 90% confidence level upper 
limit on the fractional rms amplitude of 2%. 

Finally we comment on the transition of the source 
from the banana branch to the island state. As indi- 



cated in the original classification paper (Hasinger & 



van der Klis 1989), our observation strongly suggests 
that the transition from the banana branch to the is- 
land state is most likely a continuous one (see the last 
paragraph in the last section), not an abrupt change. 
As the source moves closer to the island state, i.e., 
when the count rate decreses, the upper QPO peak 
becomes broader and its overall fractional rms ampli- 
tude becomes higher. Its centroid frequency, to the 
extent that we can tell from our data sets, follows con- 
tinuously along the expected curve extroplated from 
the banana branch. 

We would like to thank the anonymous referee for 
comments and suggestions that helped us improve 
this paper. W.Z. would like to thank Elihu Boldt 
for many stimulating discussions and encouragement 
in preparing this paper. 
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Fig. 1. — Color-color diagram of the 90 segments of 
data. The filled and open circles represent those seg- 
ments with and without obervable kHz QPO peaks in 
their FFT power spectra, respectively. The two ar- 
rows identify the two segments whose power spectra 
are plotted in detail in Figure ||. 
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Fig. 2. — Low frequency part of the power spectra 
of the three points as indicated in Figure [I]. Accord- 
ing to Hasinger and van der Klis (1989), Panels a, b, 
and c, correspond to the upper banana branch, lower 
banana branch, and the island. Therefore we con- 
clude that we have observed kHz QPOs in all three 
branches/state of the source. 



Fig. 3. — Comparison of power spectra of two seg- 
ments of data: (a) from an observation on 10 Septem- 
ber 1997 and (b) from an observation on 9 February 
1997. Although the source is in very similar spectral 
states and has very similar count rates, 3,513 cps and 
2,979 cps, respectively, they display rather different 
kHz QPO characteristics as shown by this figure. 
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Fig. 4. — QPO centroid frequencies as functions of 
count rate, where the filled circles represent the up- 
per QPO and the open circles the lower QPO. The 
two drastically different slopes are indicated by the 
hand-drawn lines. The dates on which the data were 
collected are shown. 
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Fig. 5. — Characteristics of the kHz QPOs plotted 
as functions of count rate: (a) Frequency differences 
where two simultaneous QPOs are observed; (b) The 
fractional rms amplitudes, where the filled circles rep- 
resent the upper QPO and the open circles the lower 
QPO. As in Figure |], the filled circles represent the 
upper QPO and the open circles the lower QPO. 
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